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SYNOPSIS 


A 110 M¥ imit of the Panki Thermal Power Plant, Panki, is 
selected to develop a non-linear mathematical model. The model 
incorporates all the important components of the power plant 
and is systematically presented in the text. The entire system 
is divided into three subsystems for the convenience in the 
modelling. In the present integrated boiler -turbine model dyna- 
mics of forced draft fan, fuel feeding system and furnace are 
included. Most of the previous workers in this field did not 
incorporate the fuel system dynamics. Open loop model for boiler 
unit and t'urbine unit are separately available in literature. 

In the present attempt the partial models for the boiler and tur- 
bine are also incorporated along with developed model for the 
air heating system, fuel feeding system and furnace system to 
form the integrated model. The model is simulated on a digital 
computer (DEO 1090). Extensive tests have been conducted on the 
model for various system inputs. The behaviour obtained were 
compared with those of other workers in this area. The various 
results obtained and the conclusions arrived at are presented. 

He commendations for future work are also made. 
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CHAPTER 1 


INTRODUCTION 

Tile dramatic change in the world energy situa,tion 
following the steep increase in oil prices since October, 1973 
had its impact on India also. The value of import of crude 
oil and petroleum products went up from Rs, 333 crores in 1973 
to Rs.1147 crores in 1975, creating a serious imbalance of 
payments. As an immediate step, it was decided to introduce 
measures to curb consumption of petroleum products, specifi- 
cally for non-essential uses, and identify and implement possi- 
bilities of substituting petroleum products with' other forms 
of energy, mainly coal. On the petroleum side, exploration 
activities have been intensified to locate nevr oil resources, 
particularly off-shore. In the coal sector, a comprehensive 
programme has been evolved and adopted to step-up coal pro- 
duction to meet the increased denB-nds. India is in a slightly 
better situation as far as i.:he fossil fuel resource is concerned 

In view of the fact that electricity is the most con- 
venient and versatile form of energy, the demand for it in 
India has been growing at a fast rate. Electricity plays a 
crucial role in both industrial and agricultural sectors and 
the productivity in both these sectors is closely linked with 
adequate availability of power. In view of this, power deve- 
lopment has been given high priority in the development , 
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programmes. The main sources of power in the country are 
thermal, hydro and nuclear fuels. Besides these conventional 
sources of power, the Ministry of Enerigy>^is examining the 
prospects of development of power from non- conventional 
sources such as geothermal energy, solar energy, wind energy 
and tidal energy. 

Goal Based thermal stations contribute about 58 percent of 
the electricity generation in the country. Considerable economies 
can be achieved by locating large thermal power stations near 
colliery pit heads and distributing power from them through 
extra high voltage transmission lines. It is proposed to 
initially establish, in a phased .manner, super thermal stations 
in the various regions of the country. Of the total energy, 
generated, 952,800,00 MV/Hduring 1976-77 thermal energy generated 
amounts to 572,360,00 JWH Hydro energy generated during the 
same period is 347 >910, 00 BWIjand nuclear energy generated was 
32,530,00 MWH. 

A centrally sponsored progi’ammo of establishing inter- 
state/regional links was included in the Fourth Plan and is 
being continued in the Sixth Plan also. India has now well 
connected power systems, Exclmnge of power is taking place 
regularly between a large number of states and this generally 
facilitates better utilisation of the existing capacity. 

Regional Load Despatch Centres are being established' in each 
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region to facilitate integrated operation of the constituent 
power systems. It is envisaged that, ultimately, these load- 
despatch centres would be equipped with on-line computers to 
control the power stations and power systems automatically. 

Do design efficient control systems for the power plant 
thorough understanding of the plant dynamics becomes essential. 
To study the plant dynamics actual tests on the power plant are 
highly objectionable because it is not feasible to create major- 
power system disturbances for the purpose of examining system 
behaviour and effectiveness of proposed strategies or actions. 
Simulation can be used in such a situation to judge the merits 
of these strategies under major- disturbances. Moreover sensiti- 
vity of these strategies to system changes and system model 
parameters can also be determined. Simulation provides a means 
of evaluating the system behaviour and helps in predetermining 
the critical parameters for proper design and operation and 
their responses to operator actions etc, 

])(l£iinl5?' three kinds of problems are studied in dynamic 
analysis of power systems, first kind of study is carried out 
for a few seconds of the plant operation. These are called 
short time problems. The study of instability that may occur 
when there is a fault on the genei-a ting side falls in the firs 
category. Medium time problems are of the second type. The 
study of response of the system to load changes falls under 


c*r 
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this category. The last type of problems are termed as long 
time problems. Study of start up and shut down comes under 
the third category. Present work falls under the second 
ca tag cry of medium time problem. 

The most essential ingradient to the success in the study 
of the above types of problems is an adequate mathematical 
Diodel. The model must be reasonably accurate in its represen- 
tation of the responses of interest. In the present study an 
attempt has been made to model the complete plant incorporating 
the already existing but partial models for boiler [7] and 
turbine [9J. A most general mathematical model to i*epresent a 
thermal power plant boiler-turbine unit with all the important 
subsystems are presented in Chapter 3. 

1 1 Literature Revi ew 

In the year 1958 Ghien et al.|'l] presented a mathematical 
model and its dynamic analysis for a drum type boiler. An 
oil fired, single furnace and natural circulation unit was 
considered for modelling. The model consisting of a number of 
differential equations was simulated using an electronic 
differential analyser. He conducted extensive study of the 
model and a few open loop responses were obtained. But the 
model was not capable of representing the complicated dynamics 
of the boiler. 
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A 200 M¥ uinit of the Philadelphia Electric Company, 

Crombay is coal fired, twin furnace, controlled circulation 
reheat boiler. Daniels et al. [17] (1961) developed a mathe- 
matical model for the above unit. The primary disturbance 
input was throttle valve position and the control inputs were 
coal flow rate, feed water valve position, superheater and 
reheater burner tilt positions and air flow. Unfortunately, 
the model did not produce sufficiently accurate responses to be 
suitable for control system design. In 1965, Thompson presented 
a very much refined version of the Daniels' model. The primary 
addition to the model was the use of a large number of elemental 
volumes to achieve a better approximation to the actual distri- 
buted parameter process. Although the model response was signi- 
ficantly improved, it still could not be considered accurate in 
all respects. In 1967 Eicholson et al. [16] developed an 
improved version of the model developed by Chien et al. The 
model described by Anderson ot al. [3] in the year 1968 was very 
similar to that developed by Daniels, 

IBM report [10] also presented models for drum type boiler 
and turbine using lumped parameter technique. Open and closed 
loop responses were obtained for various control parameters. 
Kwanty et al, [6] in 1970 developed an extensive model for a 
boiler-turbine system. The results of the various tests condu- 
cted on the model are presented in the report. This model is 
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relatively simple and accurate and can be used for control 
system analysis and desig-n, 

1.2 Present Work 

The entire boiler-turbine system is divided into tliree 
subsystems for convenience in modelling. Subsystem I consists 
of the, P.D, fan, air preheaters, economiser and fuel feeding 
system. Dymstpiree mathematical model for the subsystem I 
incorporates the characteristics for the P.D, fan. Empirical 
roLations have been derived to compute the fuel flow rate at 
various loads, from the actual plant data collected for a 
period of three months at Panlci site. In the furnace dynamics 
we have incorporated the effect of emissivity of the flue gas 
also. Subsystem II consists of the boiler drum, downcomer, 
waterwalls and superheaters. Turbine, rehoators, condenser and 
feed water heaters together form subsystem III, 

An integrated model for the boiler-turbine system is 
ari’ived at by coupling the tlxcee subsystems together. Models 
for subsystems II [7] and III [9] wore already existing separa- 
tely, In this thesis, model for subsystem I have been developed 
and all the models have been combined together. The entire 
system has then been simulated on a digital computer (DEC 1090). 
The behaviour of the integrated model for various test inputs 
have been compared with the results published by various 
authors. The behaviour of some important system parameters 
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such as drum pressure, flow throUi^h throttle valve, flow at the 
exit of the MP turbine, reheater steam temperature and pressure 
etc. have been obtained. These transient beluaviours are impor- 
tant in the study of sudden variations in plant inputs. 



CHAPTER 2 


DESCRIPTION OP THE SYSTEM 

2 . 1 Jntroduction 

The thermal power plant selected for the present study is 
situated at Panki. Panki is a small towiij, an industrial area, 
which is about five kilometers away from I.I.T., kanpur. The 
erection and commissioning of the two 110 MTif units have been 
done by M/s. Bharat Heavy Electricals Limited and Instrumen- 
tation Limited, Kota, 

2*2 Subsystem I 

2.2.1 Puol system 

The general lay out of the system is as given in Pig. 1, 

Raw coal from the coal mines is transported to the power plant 
by trains, Hag on tippler helps lifting and tilting of wagons 
filled with call. \Iagons are titled to 135° from the vertical 
to unload the coal by gravity. The unloaded coal is fed to 
conveyor belts, sunning over rollers, at a velocity of 2.1 m/sec. 
by electrically operated vibration feeders. After separating 
the magnetic particles from the raw coal, by huge magnets, it 
is sent through sieves. Coal particles of siae less than or equal 
to 20 mm are now fed directly into another conveyor belt lociding i 
to Raw Goal Bunkers (R.C.B.)., Gaul particles bigger than 
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20 mm sizG are crushed to the required size before sending to 
the R.O.B. There are thi'^ee R.C.B.'s for each of the 110 Wil 
units. Each R.O.B. is of 500 tons capacity. 

Crushed coal from the R.O.B. is sent to the ball mills, 
for pulveriso.tion, by Raw coal feeders. Each of the 110 FM 
unit is having three ball mills. Flue gas from the second pass 
of the furnace is tapped from below the convection superheater. 
Tapped flue gas is sent to the ball mills along with the coal. 
Hot flue gas helps pulverising the coal by heating it up. Flue 
gas is also used as the carrier of the pulverised coal from the 
ball mills to the classifier. C3,asBifier separates pulverised 
coal from the bigger particles of coal. Bigger particles of 
coal are sent back to the mills while pulverised coal particles 
are carried away by the flue gas to the cyclone separator. 
Separated pulverised coal, inside the cyclone sop^arator is fed 
to the Pulverised Goal Bunkers (P.C.B.), 

Pulverised coal from the P.C.B, is fed to the P.C. Mixer 
by the P.O. Feeding system, P.C, Feeding system consists of the 
following equipments, viz. Shut off side gate, Flap valve, 

P.O, feeder, Flexible coupling, P.I.V. Chain gear box and 
electric motor. P.C. feeder is the heart of the system. The 
body of the P.C. feedei* is of cast iron. The upper portion is 
divided into throe chambers by meajis of partition discs. The 
top chamber which is the pulverised coal receiving member has 



10 


four arms welded on a circular hub. In the middle and bottom 
chambers there are eight straight arms in each, attached to 
their respective hubs. All these three hubs are keyed to 
a Vertical slmft which rests on a thrust bearing located at the 
bottom to enable the flow of pulve3rised coal from the top and 
middle cliambers to the bottom chamber. They are so arranged 
so that the flow of pulverised coal is not affected due to the 
flow of primary air to the mixing nozzle. Though this arrango- 
ment prevents the free flow, it will help continuous flox'j- of 
pulverised coal. There are sixteen P.G. feeders which supply 
coal to P.C. Mixer at a controlled rate. The rate of flow of 
pulverised coal is varied by varying the speed of the P.C. 
feeder. The output from the feeder is proportional to the r.p.m. 
of the P.0, feeder znotor, 

2.2,2 Air system 

Atmospheric air required for combustion of the fuel is paimped 
into the combustion chamber by I'oi’ced Draft (P.D.) fan. Air from 
the P.D, fan goes to the Air preheater I and then to the Air 
preheater II. Plue gas flows inside the tubes of preheaters 
vhiile air flows outside the tubes. During the flow through the 
preheaters air gets heated mostly by convective heat transfer. 

Hot air from the Air preheater II entei'S into the Air distri- 
hution box where it gets bifurcated into the primairy air and the : 
secondary air. Secondary air directly enters into the combustion : 
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ohamber. Primary air mixes witli the pulverised coal, coming 
out of the P.O. feederSj inside the P.O. Mixer. Prom P.0, 

Mixer the mixture of primary air and pulverised coal, the fuel 
ready for combustion, enters the combustion chamber through 
burners , 

2.2,3 Firing system 

Oombustion chamber is rectangxilar in shape. There are 
sixteen ceal burners, located at the four corners of the combus- 
tion chamber tangentially along the periphery of circles, having 
their centres co-axial with the centre of the combustion chamber. 
At each corner there are four coal burners, three vapour 
burners and two oil burners. All the nine nozzle tips are tilt- 
able and are connected to a common tilting mechanism by means of 
a connecting bar. The coal nozzle tips (four numbers) are 
provided with jacket air which flows around the call fuel stream. 
The burner tips can bo tilted to fifteen degrees downwards and 
fifteen degrees upwards (i.o. a total of thirty degrees) from 
the horizontal, 

2*3 S ubsystem II 

2,3.1 Drum, Domicomer and Watorwalls [7] 

Boiler feed pump sends feed tater from the feed iciter tanl: 
into the drum. On the way from the feed water tank into the 
drum it passes through high pressure superheaters and economiser. 
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Wator collected in tlie drum flows down tlie dowicomer. Domi- 
comer is insulated and is outside the furnace, from the down- 
comer water enters the watorwall, which is inside the furnace. 
Risers, also called wateri-Kills , are arranged on all sides of 
furnace, Waterwall gets heat mostly hy radiation. Inside the 
waterwall a two loliase mixture of water and steam is formed, which 
flows into the drum tlirough headers, 

2,3.2 Superheaters [7] 

Steam is collected from the top of the drum, which is at 
saturated condition. Drum steam first goes to the ceiling 
superheater and then to the convection superheater. From the 
convection superheater steam passes through platoii superheater 
and final superheater. Steam properties are controlled between 
convection and platen superheaters, and between platen and final 
superheaters by attemporation spray. Ceiling, platen and final 
superheaters receive heat mostly by radiation, 

2 • 4 Subsystem III 

2.4,1 Turbine system [9] 

Steam from the final superheater passes through a normally 
open throttle valve, and enters a steam chest. Each steam 
chest contains four governing valves, which control the steam 
flow to the turbine. G-ovorning valves are controlled by oil 
pressure which can bo remotely controlled by a pump. Steam 
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leaving the governing valve is expanded through convergent 
divergent nozzles and attains very high velocity. High velocity 
steam gets expanded in two rows of moving blades called the 
governing stage. The partially expanded steoin then gets expanded 
in High Pressure (H.P.) turbine. Steam from the H.P. turbine 
is sent for reheating into Triflex rehoater followed by exit 
reheater. Reheated steam is then expanded in Medium Pressure 
(M.P.) turbine followed by Low Pressure (L.P.) turbine. The steam 
exiting from the l.P. turbine is condensed in two condensers 
each,, and the condensate is passed tlu^ough a series of L.P. feed 
water heaters. 

2.4.2 Feed water heaters and economiser [9] 

Condensate, feed water from the condenser is sent tlirough a 
series of five L.P, feed water heaters. Inside the L.P, feed 
water heatei' condensate gets heated by extracted steam from L.P. 
and M.P. turbines. Hot oondonsate at the end of L.P. food water 
heater enters the Feed water tank. Boiler foed pump pumps the 
food water from the tank to H.P. foed water heaters I and II, 

Feed water gets heated inside the H.P. hecLters by the extracted 
steam from the M.P. turbine and by a small amount of exit steam 
from the H.P. turbine. Foed water from the heater now enters the 
economiser. Inside the economiser feed water gets heated finally 
by fluo gas. Hot feed water then passes to the drum. Further 
details of the system are given in Appendix I, 



GILiPTER 3 


mTHEmTIGAL MODELS 


5 • 1 Int roduotion to Modelli ng 

In an era of inoreasing electric power utilizations, plant 
availability has becomo of paramount importance and consequently, 
its ability to ride smoothly tiu''ough major as well as minor 
system disturbances has become essential. Porliaps the most 
essential stop towards better understanding of the gonerf^tion 
system and achieving the above is obtaining an adequate matho- 
rnatical model. Roughly speaking, the basic roquireraents of 
such a model are : 1) it must be reasonably accurate in its 
representation of the responses of interest, 2) it must bo 
relatively simple, and 3) it must bear close resomblencc to the 
physical process, i.e., model parameters should be 
readily related to physical parameters, 

¥e Inive adopted a theoretical approach in developing tlio 
system rnodGl, In this method development of the model is based 
on some physical laws and also on some seni-ompirical relation- 
ships. Yariables like temperature, velocity etc. that depict the 
behaviour of components, like waterwalls and super heaters, are 
continuous functions, of tiiae and spaco. Dynamics of such 
components am be adequately described by partial differential- 
equations , Siiiailation of such a model nocesaitaies considerable 
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computer time which becomes prohibitive in view of the large 
number of coraponents to bo inc or pore, ted in such a model. The 
entire system is divided or lumped into several sub-systems, 
across which changes in physical properties take place. The 
clnngos across aich of the blocks are represented in the form 
of nathematical equations b3^ applying momentum, laass and onerg^?' 
balance. The mathematical model developed is presented, systo- 
inatically, in the following articles, 

3*2 Model for Sub sy stem I 

Subsystem I consists of economiser, f.D. fan, air pz-eheater 
I, air preheat or II, fuel feeding system, and combustion chamber. 
Pood water pumped out from the food I'/ater tank passes through the 
high pressure heaters I and II and then through the eoonomiser 
before it enters the boiler drum. Atmospheric air utilised for 
combustion of fuelling is sent into the combustion chamber by tho 
P.D, fan. Air coming out of the P.D. fan passes through air pre- 
heaters I and II where it gets heated by flue gases as explained 
in Chapter 2, Tho block diagram showing tho flow path is shown 
in Pig. 1. 

3.2.1 Economiser (watoi’ side) 

The following assumptions were made ; 

1, Average values of the properties of the water passing 

through tho economiser are astsumed to be the exit values 
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2, Water entering the econoniser from the high pressirre heater 
II has the instantaneous properties of the v/ater leavjjig 
the high pressure heater II, 

3. Water passing throiigh the economiser reoeives heat mostly 
by convection. 


The rate of heat transfer from the metal tubes to the viator 
is calcuhated using omperical relations. The flow of -v/ator 
through the economiser is highly turbulent, and the heat transfer 
relations are of the type [2,3] 




" \m2 


(3.1) 


where 


^EI'W 

%I2 

^EM 

^E¥ 


= heat transfer rate from the oconomiser tube metal 
to the water (kcals/soc. ) , 

= a constant to be determined at the tine of 
initialisation, 

= feed water flovr rate through the oconomiser 
(kg/soc.), 

= temperature of the oconomiser tube metal ( c), 

= temperature of the feed i/jater leaving the 
economiser ( cj. 


Initial value for can be calculated by the following 
energy balance equation [4] utilising the data for the resident 
heat content collected directly from the plant. 
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SlWl “ ^^B¥ “ ^HH2^ 

whore 

^BV/ ~ specific heat of econoEiiser -vvator (kcals/kg ^c), 

= resident mass of v/ater inside the economiser (hg), 
~ temperaturo of v/ator coming out of high pressure 
heater II (°c). 


3.2,2 Economiser (gas side) 

Sihiilar assumptions (as in 3.2.1) hold good on the gas side 
of the oconomisor. 

Energy balance on tho motsil tubes of the ccoiiomiser lead to 
tho following equation. 


Q. 


BGM 


Q 


BM 


= 0 . 


BM 






El'h 


(3.3) 


where 



■^BM 

^BM 



= heat transfer I'ate from flue gases to the 
oconomisor metal x-rall (kcals/sec. ) , 

= temperature of the oconoraiser 'bubo wall ( c), 
= specific heat of the economiser tube raotal 
(kca.Ts/kg °c), 

= mass of tho economiser notal tubes (kg). 


Hoat transfer rate from flue gas to tho economiser metal 
tube can be ccxlculated using tho relation [2,3]? 



0 ^ 
^BGM 





(3.4) 



whore 
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^ constant to be dotorrained at the tiiac of 
in it ia 1 i s a t i on , 

= flow rate of flue gas through the ecoiioinisor 
(kg/sec. ) , 

^EG- " tenporaturo of the ccononiser fluo gas ('^c). 

3.2.3 E.D. fan 

Specific assumptions made during; tho raodolling of F.I), fan 
are tho following ; 

1, Volume of air flowing out of F.D. fan is linearly related 
to the external load on the power plant. 

2, Atmospheric conditions exist at the inlet to tho E.D, fan. 

Tho follovring ompirioal formula has boon derived from tho 
plant data for calculating the rate of air flow across the E.D, 
fan , 

V^p = 2818.0 X LOAD + 110.0 x 10^ (3.5) 

whore 

V^.^p = Volume flov/ rate of air across tho F.D. fan 
(n^/hr) , 

LOAD = exteriml load on tho powei'* ploait (MI'f). 

Mass flow rate of air a, cross tho F.D. fan is given by 

^’^AF " flF '^AF 


(3.6) 



where 
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¥ 


AF 


(j 

af 


Y 


AF 


Mass flow rate of air across the F.D. fan 
(. kg/sec. ) , 

density of air coning out of F.D. fan (kg/cr), 
volume flow rate of air across the F.D, fan 
(m-^/sec. ) , 


3.2.4 Air preheater I (air side) 

Similar assumptions as in 3.2.1 hold good in the oase of 
air preheater I also. 

Heat transfer rate from the tube notal x^all of the air 
preheater I to the air can bo foimd out using the formula 


Sim “ Sima ^'^af^ 


PIM 


SlA^ 


( 3 . 7 ) 


where 


Q. 


'Pim 


SlM 

^PIA 


G 


PIFIA 


Heat transfer rate from the tube r.ietal to the 
air (kcals/sec. ) , 

temperature of tiio iJX’cjhoator I tube metal (°c), 
temperature of the preheater I air (°o), 
a constant to bo clotorminod at the tirae of 
initialisation. 


Initial value for Pp2.A calculated by the folloxving 

energy balance equation [4] u.tilising tho data for the resident 


heat content Qp^jj^i colle°'t'2ii directly from the p.lant. 
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SlA ^'^PIA “ ^AP^ ^%1A 


(3.8) 


where 


M. 


G -, A 

plA 

'^AP 

■PIA 


Speoifio heat of the air passini^j through tho 
preheater I (kcals/kg ^c), 

temperature of the air coning out of the P.D, fan 

( c) , 

resident mass of the air inside the preheater I (leg) 


3.2.5 Air preheater I (gas side) 

On x^riting down the energy balance for the tube metal of the 
preheater I wo arrive at tho folloxfing ecjuation 


Q-, 


PIGM 


SlM “ SlM ^%1M dt 


• ( T ) 


(3.9) 


where 


*^P1GM “ Heat transfer rate from flue gas to the preheater I 
metal (Icaals/scc, ) , 


PIM = specific heat of the preheat or I tube metal 
(kcals/kg °c), 

mass of the preheater I tube metal (leg), 


M. 


PIM 


^PIM ~ temperature of the preheater I tube metal (°o). 

Heat transfer xmto from tho fluo gas to the preheater I 
tube metal can be calculated using the formula 


PIGM 


'^PIGM 


PIG ^PIM^ 


( 3 . 10 ) 
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*^P1GM ~ ^ Goiistant to loo detorrained at tlio time of 

init ia livsat i on , 

fpiQ. = tenporaturo of the flue gas pavssing tlnrough the 
preheater I (°o), 

= mass flox^ rate of flue gas passing across the 
pr ohea ter I ( leg /sec.). 

3.2,6 Air prohoater II (a.ir side) 

It is assumed that the air leaving from the air preheater I 
and the air entering the air preheater II has got the same 
properties , 


Heat transfer rate from the air prohoater II tube metal to 
the air passing through the prohoater II is oilculatod using the 
following erapirical relation. 


0 , 


P2MA 


n ^ v 

^P2Mil ^-^121 


^P2il^ 


(3.11) 


whore 


%2mi 

^P2MA 


T 


P2M 


I 

^P2A 



= Heat transfer rate from preheater II meta3. to 
air (kcals/sec, ) , 

a constant to ho determined at the time of 
initialisation, 

= temp el's turo of the prohoater II metal ("^c), 

= teraperature of the preheater II air (°g), 

= rate of air flow through the preheater II 
(kg/seo. ) . 



Initial value for calculated by the following 

energy balance equation [4] utilizing the doata for the resident 
heat content oollectod directly froi.i the plant. 


Q. 


'P2MA.1 


M. 


G. 


■P2A ^P2/l ^-"P2/i 


('A 




(3.12) 


where 

Mp 2 A = resident nciss of air inside the preheater II (hg), 
^P2A " specific heat of air inside the preheatex' II 
(Iccals/kg *^c), 

IpiA = tempera.ture of air leaving preheater I (°c). 

3.2.7 Air preheater II (gas side) 


Energy belonice on the prohoater II tube notal le.ads to the 
following G qua t i on . 


^P2GM " SeMA 


^P2M ^P2M ^'^P2I'P 


(3.13) 


whore 


S2GM ~ 


^P2M 

M 

^T2M 

^P2M 


heat transfer rate fron the preheater tube gas 
to the raotal (kco.la/soc, ) , 

specific heat of the prehuutor II tube metal 
(k cals /kg °c)j, 

mass of the preheater II tube metal (kg), 


tenpei’ature of the proheater II tube metal 
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Heat transfer rate from flue gas to the air preheat or II 
tube metal is calculated u.sing the orapirical formula 

S2GM ^^P2G ~ '^'P2M^ (3.14) 

where 

^P2GM ~ constoait to be detoriainod at the time of 
initialisation, 

Tp2Q = tempera turo of the preheater II gas (°c), 

3.2.8 Fuel feeding system 

Functioning of the fuel feeding system is crrplainod in 
Chapter 2, To derive an empirical relationslifp between the 
external load on the power plant and number of pulverised coal 
feeders (P.G, feeders) which are live for that particular load 
wo resorted to the available plant data. 

Readings of load against tine have boon noted continuously 
for throo months. For tho above period the nuubor of live P.O. 
feeders were also noted against tirae. All the readings have boon 
suitably tabulated and ultimately a graph showing the number of 
open P.O. feeders versus load in IPd have been prepared as shown 
in Fig. 2. 

The following empirical rolationship is derived from the 
graph between load and the number of open P.G, feeders. 



PqP == LOAD/21. 0 + 9.095 
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(3.15) 

where 

PqP = No. of live or active P.O. foedors, 

LOAD = External load on tho power plant (Mh). 


To find out the fuel flow rate into tlivi conbustion 
chaiaber the following eq.uation is used. 


- ■^PGF ^CF 


(3.16) 


whore 

V/j, = Fuel flow rate into tho combustion chamber 
(kg/soc. ) s 

WpQp = fuel flow rate through one P.O. foedor (Icg/sec.), 
3.2.9 Combust ion clmnibor 


Tho combustion chamber is asswned to bo a black onolosuro 
in xirhich radiation exchange between the entire gas volume and 
enclosure boundary takes place, Fioan beam length approach [5] 
is used to analyse the situation, Fieaii boa:;, length is defined 
as the required radius of a gas honisphero such that it radiates 
a flux to the centre of its base eqoial to the average flux 
radiated to tho area of interest by tho actual volume of gas. 

Flean beam length can be used as a characteristic dimension of 
the gas volume and regarded as a constant. The recommonded 
equation to calculate the mean beam length is given below [5] 
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= mean beam Icaigth (m) , 

^PUE ~ effootive furnace volume (ix), 
j'^pUR “ effective furnace area (m ). 

Eniissivity charts have been developed by Hottel [5] from 
many experimental nioasuroments . The flue gas may contain 
car bond i- oxide j x^ator vapour carbon liionoxido, methane and a 
few other gases. But it is assumed that oxcopt water vapour 
and oarbondi-oxide, all' are spectrally transpon-ont gases. It 
is also assumed that the gas is ]an.ving uniform temperature insid 
the chamber. 


0-.r 


Id 


Graphs aro available [5] for q vs IC and Eqq 

t.tI'I /%v» rT\ O' o vn c* 

COg 


vs T_, K, xdiore Sp. ^ and e^,r^ are the omissivities of water 


gw 


“H2O 


vo.pour and curb ondi- oxide respoctivoly . T„„ is tho gas tem- 

perature insido tho chamber. These give the t'otal emittance 
of carbondi" oxide and water vapoixr obtoinod exp or imen tally 
using a mixture x^rith non-absorbing gases so that total pressure 


of the gas was 1 atmosphere xdiilo partial pressures of carbondl 
oxide and water wore varied. 


Per a gas mixture other than 1 atmosph 
broadening correction is to be applied [5]. 
a multiplying constant ^002 case of 


oro, a pressure 
This is given as 
ca r b end i- ox id e . In 
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tliG case of water vapour the eiaittanco is influonccd in a 
slightly nore Gomplex mnnor by both the rjartial pressure of 
the water vapour and the total pressure of the gas mixture. The 
correction for water vapour q is dependant on the moan boara 


length also. 


for our problon ho. and Gjj q liave been estimated from the 

graphs [5] vs t P^^q + Pqq /2.0, Cj.j. . vs 

tota l ^H 20 ^ i l 


"2 ^ ^H20 ^'c» 


where 


Ptotoi ^ total pressure of the gas, 


P 


H2O 

'00^ 


= partial pressure of water vapour, 

= partial pressure of carbondi- oxide 


^ 00 , 


and Og. Q can be assumed to be constants in our case for the 


partial pressures inside tho combustion chamber are assumed to be 
constants . 


As in the present case if COg and water vajjour are both 
present in the gas mixture, an additional quantity Z.\ e must be 
included to account for tan emittanco reduction resulting from 
spectral overlap of tho carbondi- oxide and water absorption 
bands. This correction factor is found from the graph [5] 


b.^z vs q/(P 


00, 




for a mixture of carbondi-oxide 


and water vapour in a non-absorbing carrier gas the eraittance 


is then given by 
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- ^00, 


'GO2 ■*' ^H20 %2*^ " 


(3.18) 


In tiiG prc'Sent case oiipirical r olatioiishipc have hcon 

derived between and T as well as ,, and I . The 

eU 2 gw ri 2 G gw 

equations are presented below. 


= 0.29 - 0.0001 + 273.0) (3.19) 

= 0.212 - 0.05 (Tg^.^ + 273.0)/600.0 (3.20) 

where 

1™,. = gas temp)oraturo inside the chamber, 

gw 

The process of combustion of pulverised coo^l in a furnace 
is extroiiiely complex. A macroscopic view of the process indi“ 
cates that if the rate of flow (kgs. /sec.) of coal having 

a heating value equal to is burnt, the gas teraporature with no 
heat loss [6] would bo 


T 

gw 


Wp Op pp 


'■^P2A 


( 3 . 21 ) 


where 

pp = efficiency of burning, 

= specific heat of the gas, 
b* 

= mass rate of flow of flue gas through the furnace 
(kg/sec,). 

Writing dow tho mass bala.nco inside the combustion chambor 
we get 
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Wg = + Wp (1 - (3.22) 

where 

^'^ASH ~ content in the pulverised ccc.ly 

''^AF ” Class flow rate of air into the conhustion chanbor 

(kg/s GO. ) . 


The ceiling and convoetion superheaters arc clubbed 
together and is torned as primary superheater. This forms the 
first section. Platen and exit superheaters togetlior is also 
considered as a single unit called the secondary .suporhoator. 

Both waterwalls and secondary superheater lie in the first 
pass of the furnace while the primary suporhoator is in the 
second pass of the furnace, which rccoivos hoat energy mostly 
by convective hoat transfer. Secondary superheater and water- 
walls receive heat mostly by thermal radiation from the flue gas. 


Taking the heat balance on the gas shows that the not 
average heat flux being I'emovcd at the w:ill i:-j the gas emission 
minus tho amission from the toII that is absorbed by the gas. 




( 3 . 23 ) 


where 




= heat transfer rate from fluo gas to tho viatorwalls 
(kcals/sec. ) , 

2 

= effective watericrall surface area (m ) , 
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’^¥11 

= waterwall metal tempera tux'o (°c)j 

= absorptanco of the gas for radiation 

emitted 

a 

from the wall at temp era, ture Tyjyj? 

= Stefan-Boltzmann constant (kcals/sec 

2 0, 

m k ) . 


Similarly for tiio secondary suxJeriioat or , tho radiative 
heat transfer rate fi*OLi the flue gas to tiio metal walls can be 
oalculatod by the followi]ag equation 


Q 


SM2 


= "^32 <'^g^V 273)-^ - a^. 


g 


273 )"^] 


(3.24) 


where 


*"^SM2 “ heat transfer rate from flue gas to the secondary 

superheater metal (kcals/soc. ) , 

Ag 2 - effective surface area of the secondary superheater 
raotal exposed to the furnace (n ), 

^^SM2 ~ ■tc^rnporaturo of the metal tube of secondary supor- 


heater (°c). 


Duo to tho considerable boat transfer rates into the mtor- 
walls and primary superheater in tho first pass of the furnace 
tho temperature of the flue g-as docroasos as it leaves the 


first p.ass. 


An energy balance on flue gas in the first pass 


is done to find out the temperature of tho flue gas .leaving the 
first pass. The derived equation is as follows. 


■G¥P 


[(Tg^-i. 273)" 


(Q 


. “COBP 


) 1/4 


273 (3.25) 
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^Ct¥P ” ■temporaturo of the fluo ga-3 leaving the first 
pass of the furnace (°c)j 

GOEP = a constant to bo doterminod at the tine of 
initialisation. 

Heat transfer rate fron the flue gas to tiio primary super 


iioator can be oalculatocT using the empirical relation, 






0.6 

G 






vriiorc 

^i-ms 


'P 

■^spa 


= heat transfer rato from the flue gas to the 
pr imarjr superheat or ( kcals/ s oc . ) , 

= a constant to be deterrainod at the time of 
initialisation, 

= tomperaturo of tho primary siiporheautor nots.l 

( c). 


5 • 3 Mod el for Subsystem I I 

Subsystem II consists of tho drum, downcomer, wa,torwalls , 
primary superheater and secondary superheater [7]* All tho 
variables are defined in the Nomenclature. 

3.3.1 Di'um-Downcornor-'watex’wrills 

Drum-bowncoHor-watei’wallo system is considered as a 
single loop in which fluid f low , is due to tho density cixangos 
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only and henco natural circula-tion enists in tlio loop. This 
systoD is convoniontly divided into tliroo separate sections 
viz. drum, downcomer and watorwalls to facilitate lumped 
paraootor reprosentation [7]. ' 

3.5.2 Drum 


Sp G G if i c ass uiiip t ions s 


1. The drum is well insulated and heat energy stored in the 
drum mterial is negligihlc, 

2. Outlet flow rate of stoan from the drum is assmed to be 
proportional to the square-root of the product of drum pressure 
and density [8]. This is made possible by assuming an imaginary 
tlirottle valvo .at the outlet of the drum. 

5. Food wator is supplied at saturation tomporature to the drum 
and variations of propertios of vjator ?.nd steam within the drum 
arc very small. 


imss balance and energy balance on the drum water yield 
tho following equations ; 


^^‘^0 0 *“ “ri 


and 


^'Mi^ 9-‘&' ^ drw ^drw ■' drs ^drs 


(3.27) 


^Vi \'i 




^"eo + "Vo 

- V h -I- d V h ) 

~ dT ^ drw ^drw "^drw drs drs Vrs-' 


(3.28) 
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Following empirical relations are used to deternino certain 


variables 

V, = 
dr 

which appear in the above equations, 

Y 4 - V 

drs drw 

(3.29) 

drw “ 

^4 ^5^^dr ^3^ 

(3.30) 

— 

drs 

^12 ^13 ^^dr ^3^ 

(3.51) 

^Vrw “ 

^1 + ^2 ^^dr ^3^ 

(3.32) 

^^drs ~ 

^10 ^11 ^^dr ^3^ 

(3.33) 

whero Og, 

, Oil and 0^^ are constants to bo del 

;erminod by 


using linear fit for the enthalpies and densities from the 
steam tables. 

The stovani flow rate out of the drum is ailculo.ted by using 
the formula 


^‘Vi *^throttlo ^^dr ^drs^’ 


(3.34) 


whore 0,, = a constant determined at the time of 

thrortle 


ini t ial isa t i on , 


3.3.3 h own c om or 


Mass, energy and raomontum balance on the circulating, hot 
water through the downoomer yield the follov^ing eiuatioiio : 




(3.35) 
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^■^do ^do^ dt ^^d ’ 


(3.36) 


and 

r di 

_ I : 


-f T 2 2 

P ) - 11§L . i “ 

do g ^1“ 2g 5 


di 

di 


s^d ni :. 


do 

2-75 — 
d “’do 


■] “ -^di gir h 


^d 


(3.37) 


Empirical relation used to oaloulate the density of 
down com or outlet wator is given by 


?do = “le^^do + + ®17 (‘^do - Us) + «19 


(3.58) 


3.3.4 Watorwalls (wator side) 


The mass, energy and momentmi balance for the water in the 
waterwalls are as given boloxir : 





m. 


wo 


^wo^ (it' ^^w \fO^» 

Wo + W) = sir (W ^’wo Wo)' 


and 

- Wr) = 


ra 

'[ — ‘2 


2 

do 

■“^do 



2f L .p 
W ^ JJX 


2 

1 “wo 

D 

W w 

p ^ 

Jwo 

^ ? 

. 2 

+( 

j wo 

+ 1 
©i 

, , 

1 

9 

-* W0>, T- 

2 


+ 



d_ 

d^ 





(3.39) 

(3.40) 


(3.41) 
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Enpirical relations used to oilculate the two-phase 

friction factor, f^, density of the two-plmso mixture leaving 

the waterwall, *5'^^ and the entlialpy of tlw two-phase mixture 

leaving the waterwall h ^ are given helovm 

wo 


L'l + 2400 


(3.42) 


j \ro 


p 


+ 


‘ wot 




?wow 


(3.43) 


and 


h 


wo 


X h 


wos 


+ (1'- X) h 


wow 


(3.44) 


3.3.5 Wa t oi’wa 11 s ( ga s side) 

Energy balance for the wo.terwn,!! tube metal yield the 
following aquation .♦ 


~ ^'WM 


(3.45) 


Heat transfer rate to the two-phase mixture from the 
watoi-wall is calculated using the equation 


= CONST (T^,^jvi 


^dr ^ 


(3.46) 


An empirical equation is used to calculato the tomperetturo 
of the liquid in the waterwcills, which is given by 


T 


dr 


14 + I 5 (Ir + 1) 


( 3 . 47 ) 
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3.3.6 Superheaters 

Superheaters are represented bj>- throe sections. The 
coiling and convection Bujperhoators tegethor is torned as 
priiuary superheater, which is the first section. The third 
section is the combination of platen and exit suporhoaters and 
is called the secondary superheater. The second section is 
assumed to bo a unit which affectively simulates the flow 
dynamics of the superheater lines. This c:in be imagined as a 
connecting link between the primry and secondary superheaters, 

3 . 3.7 Primary superheater (steam side) 

Eiss and energy bOwlanco relations for the steam passing 
through the primary superheater is exprossod as 




^’^0 ^'^DSl 


) = 


h '-'si 




(3.48) 


and 


d_ 

dt 


O-igi) = L'«si + ‘-Vi ‘Vo£ 


1 


^■’^0 ^^sl '^DSl ^siT^sl 


(3.49) 


The oxi^rossion for clrxnge in enthalpy of sujjorhoatod steam in 
the primary superheater can be written as 



d / c; 'v -I 1 


'sl^J ■(TTlk^j)' 


( 3 . 50 ) 


The pressure in the primry superheater is calculat-ed -using 


the GCLuation 



q1 


^^ 2 ^ '->.31 


56 

(3.51) 


3.3.8 Primary sup or heater (gas side) 


An Giiorgy baloaicG for the priraarj/' sup erh carl: or yield the 
following equation ; 


(*^siya “■ ^“Si^ “ ^si ^^1 he ^^SHi^ 


(3.52) 


The rato of heat transfer from tubes to the stoan is 
calculatod using tho expression 


Qo 


SI 


(I 


sm 


^sl^ 


(3.52) 


where is caloulatcd using tho equation given below 


)1 




afl' 


"Rl^ ->sl 


^ sll 


) -1- Tc 


T1 


(3,53) 


3.3.9 Connect ing link between prinrvry and .secondary superheaters 

It is assumed that in the steady state flow is proportional 
to tho pressure drop between the superheaters. Thus, we can 
write : 


T 


GL dT 


(Wq) 


P 


(Pc, 


Si 


s2' 


V/, 


(3.54) 


whore Tqj^ 


tine delay constant for the- connocting link (soc. ) 


3.3.10 Secondary superheater (steam side) 


Mriss and energy balcmce equations for the steara flowing 
through the secondary superheater are, presented below. 
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1 ''DS2^/''b2 = h 


— W-| + W-p.qo)/V 


( .)' s2^ 


(;>.55) 


and 


d 


dt^^s 2 ^ - LQ 32 + ''•'^■n.qo 


'0 ■’■ ^'DS2 “ ^^1 ^s2l 


'1 

is2 ^s2 


(5.56) 


rc'Ssion for the rate of cnang'o of enthalpy of steonn ii 


by 

d 
d 


dn:'^^s2^ “ -''' 32^3 rT'^TcTjT (5.57) 


J P 


s2 d 

5^ - s2^j ri'^Th' 

'S2 b 

Tho pressure 1^2 bs oal epilated fror.i tiio reJ-ation 


^’ s 2 “ ^^^1 ^ s 2 “ ^^2 3 j ' s 2 


(3.58) 


3.3.11 Secondary superheat or (gas side) 

An onorgy balance on tho secondary suporlioater metal yield 
the following equation ; 


(Q 


Qqo) ~ 8,qo Mo (l.qMo) 


SM2 


'S2 


S2 “2 dt ^-^SM2' 


(3.59) 


o rate of hooit transfer from the tu.bos to tho stean is 
calculated using tho er.ipirical relation 


"“S2 


"HS ^"1 


(¥-, )'^'® ('Sc5ni2 “ ^s2^ 


(3.60) 


for tho superheated stoara the temp ora ture is calculated 
using the expression [7] 


'■‘■S2 - %2^^^s2 ” 


sl2 


^ %2^ -^2 “ ^’sl2^ ■'* ^ 


312 


(3.61) 
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Tho flow rato of steam from superheater tlirough the 
noriially open throttle valve is calculated using the expression 




1 




(3.62) 


3 • 4 Model for Subsystem III 

Subsystem III consists of the governing stage of the 
turbine, high pressure turbine, reheat ors, nediu-n pressure 
turbine, low pressure turbine, feed water heaters and condenser 
[9]. The model is presented below, 

3,4.1 Governing stage of turbine 

Thoro are four throttle valves ['9] to control the admission 
of steam to the high pressure turbine. These four valves are 
followed by a nozzle ring followed by two rows of moving 
blados with an intervening row of fixed guide blades. Those 
coriiStitute the governing stage of the turbine. Two of those 
four valves are adjusted to open at the sane oil pressures and 
so both the valves are lumped together into one unit. As a 
X’GSult the steam admission Ins been nodellod through throe 
valves instead of the four vo.lvos. The assumptions made at 
each stage is given in [9], The mass flow rate through the 
throttle valve is given by 

= -Q™-- Jgg (pyg “ (3.63) 
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Mass flow rate thxougli tho riozzlo ring is given by 


(K) 




G f 


THI ^^THI 




(3.64) 


Assuning tlmt there is no storage in the nozzle ring the 
flows given by eqns. (3.63) and (3.64) are equal. Hence for a 
particular load we can derive the expression for tho ra.to of 
flow of steam as 


(W,_j_J' 


A ^ G V 


a2 e 

A rp J gg 


P 

'T -*11-11 




?gr. (P 3 C 
DO OD 


Pjjp) 


(3,65) 


If the pressure ratio across the valve and nozzle con- 
bination is loss than the critical pressure less 

than 0 , 5054 ) then senic conditions exist and tho flow is inde- 
pendent of tho down stream conditions. In that case flov^ 
equation can bo expressed as [6] 

Pc 


Wrji “ Aqi 


■ss 


Vloc; 

D D 


( 3 . 66 ) 


'Iho sta'bo rela-fcions used for do-bermining nozzle ring steam 
temperatures obtained by using empirical fits from stoo-n tables 


arc given below, 

T^Hi = + 1.g3 

\'H2 ^TH2 ^%2 

'^TH3 ~ ^TH3 ^%H3 


620.0, 

(3.67) 

620,0 

(3.68) 

620.0 

(3.69) 


and 
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are 


ToiipGraturoa of G-teau streans 
calculated using the following 

Pt 


n-1 


T 


HPl 


1 


THl 


'i "D y ) 


THl 


P lirl 

m _ m / HP N n 
n-IP2 “ •^TH2 


at the exit of the noaalo rings 
equations , 


(5.70) 

(3.71) 


and 


1 — ip 

•H-IP3 ■'■TH3 


P Sri 
^ TH3 


(3.72) 


Following Gi'ipirical relations are used to detoruino the 
enthalpj.es of the exiting stean fron the nozzle ring 


%l?l ^E£ ^HPl 527.6, 


(3.73) 


%IP2 = ^PIP (3,74) 

and 

Hj.jp 3 = -0.3 Ppjj, + 0.617 + 527.6 (3.75) 

3.4.2 H igh press ur o t ui* b in e 

High pressure turbine irapulse blades are assuraed to 
ideally convert the kinetic enorgy into work. Further, it is 
also assumed tPiat the relationship between enthalpy, specific 
vo3.urao, and pressure for the superheated stean follovrs Callender's 
empirical relationship for all non-constant prossuro processes 
[ 10 ]. 
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Jho tlo¥ througii tho liigli pressure turbine is related 
bo the pressure drop across it. The expression used is 

2 


.V/,) 


^^1 HP ^^HP '* ^ 


OR 


) 


(3.76) 


The rate of change of outlet floxj of H.P, turbine with a 
change in tho inlet flow into the turbine is expressed as 


d 


d^:' 


(1 


^o) "4^ 


W 


OR 


(3.77) 


Expansion of steam in the turbine is assxmod to be isentropic 
(n = 1.3). Temperature and ontlialpy of tho exiting steam from tho 


H.P. turbine are given by 

n-1 


Tn 


OR 


T 


and 


- 0.5934 Pqj^ -t- 0.65 + 534.069 


(3.78) 


(3.79) 


The overall, offioiency, 'Hjjpi- of the H.P, turbine is 
corrected by a faotorj rip, to got the actual off icienoy, 

’’^HPA' turbine. Correction is done to conponsatc for the 

effects of stage leakages, root and tip interference losses and 
rotational losses which increase with the flow rate. Actual effiei 
offieicncy and the ;;\cjtual reheat enthalpy are then given by 


^IIPA " ^HP ^P 


( 3 . 80 ) 



42 


where 

and 


ri 




- w 

TWoj^A/^Crrw 


(3,80A) 


“ %IP “ ^HPA " %-i^ (3.81) 

The power produced "by tho H.P. turbine in terns of tho 
flow rate and enthalpy drop in the H.P, tun-bine is expressed 


as 


POW^ = 0.0042 


(E 


HP 




(3.82) 


3.4.4 Triflex roheator (stoan sido) 

Itiss, nonentum and energy balance on the stoan side load 
to tho following oquationo ; 


¥ 


0.R 


¥. 


ROl 


_ cl ( i-, 
" ( ■? 

2 

wf 


Rl ^R1 


), 


OR 


^ROl 


OR ^'1 

"ri r- 


V/n-, 


OR 


(3.83) 

(3.84) 


and. 


Y h ^\l \01^ = ^'Wl “ %L01 ^'^ROl '^OR ^5. 85) 

Tonporaturo of the stoan in tho Trixlox reheator and the 
heat added to tho steam from tho metal walls are given by 

'■^RSl " °*^®5 (Pjyo1 “ 24.47) + 1.81(I-Ijjq3^ - 772.0) 4- 397.0 

(3.86) 


and 


Q 


RSI 


‘^RMSl 


, 0.8 


ROl' 


(1 


RMl 




(3.87) 
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A large volume of steam is occupied “by tlio rehoater a.nd 
so there is a storage effect in it. As a result we can. con- 
sider it as a tank across which, the changes in flow are 
approximated by a first order transfer function. 


The following equation is derived to represent the rate 
of change of steam flow leaving the Triflex rohoater. 


d V _ '"^GR J ■'^RI2 

clT ■“T';^ 

3,4.5 Triflex rehoater (gas s 

The energy storage in the 


( 3 . 88 ) 

idc) 

Triflex reheater leads to the 


f oil owiiig equa t i on . 


Q, 


HGl 


Q 


“RSI 


= 0 


TMl \ 


d 

df 


^%M1^ 


3.4.6 B.xit r cheat or (steam side) 


(3.89) 


Mass, laomontum and energy ba,lance equations for the 
steam in the exit roheator are 


\l2 




d 


.R02 = au 


^ 4r2 ''^R2)» 


\l2 . ^2 d 


P-rn-o P-nno - I'DO — ^^^RI 2 ^ 


■RI2 

and 


■R02 " ■"R2 ,...2 

•R 


“2*^ 


(3.90) 

(3.91) 


1 d__ 
Y dt 


%2 % 02 ^ 


°“RS2 %I2 %I2 


\02 ^’^R02 


(3.92) 
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llic empirical relations used for detormining exit relieater 
bciaporatur 0 ana bho boat added to stcan in tlio exit roiieator are 

■^.RS2 ~ ^^R02 ^ 23.12) -f- 2.05 (^■Ij^02 “ + 535.0 (3.93) 


and 


Q 


“RS2 


^RMS2 ^%02^*^*^ *^^R112 “ 


■RtJ2' 


(3.94) 


AS ill tliG Ci'iSG of Triflex rGiicator, the rate of cliango of 
f].ow of stoain out of the exit reheat or is giv on by 






^’^RI2 " ^'^R02 


T 


8 


(3.95) 


3.4.7 Exit relieater (gas side) 


The energy storage in the metal tube mil of exit rohoater 
is given by the following equation 


S.G2 " 


^TM2 ^2 ft ^'^RM2^ 


(3.96) 


3,4.6 Medium pressure turbine 


Inlet flow to the M.P. turbine is supplenonted by attain- 
poration flow to oontrol the tempera turo. This is o,lso incor- 
porated in the model. Inlet flow into the I'l.P. turbine is given 
by 


“ ^^R02 ^ '’^ATMP 


(3.97) 


Energy balanco for atteiiiperator is given by iih..' exprossion 


'102 l02 = “mi '^mi ''Aim “am® 


(3.98) 
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Tho flow through each turbine is related to their res- 
poctiVG inlet and outlet pressures by tho following expressions. 




■R2 ~ 

(3.99) 

2 

w 

•qvipi 

= ^3 


(3.100) 

"^rp2 

= ''4 

i ]v£p2 ^^]yD?2 “ ^riP3^' 

(3.101) 

'^MP3 

= Kfr 

5 

•'lCP3 “ ^MPO^* 

(3.102) 


The rate of cirnigo of flow in each of tho turbine inlet 
is described by, 

h(W,^rpT ) = 




■(^IVTOo) 


d 

dt^”MP2 

d 




dt'**]yiP5 


and 


('^'^TVTPn^ ~ 


(1 - 

^4) “ 

Wj^l 

( 3 . 103 ) 


^2 

■■ f 

(1 - 

O5) “ 

Wj,jp2 

( 3 . 104 ) 




(1 

^6^ ^MP2 


( 3 . 105 ) 



(1 - 

O7) 

m 


(3,106) 


dt^"MP0 

J 

Tomporatures at each stage of tho turbine are determined 
using the following equations 

T. n-1 




MPl -"tlPI'P 


m nJgl) II 


(3.107) 


MPI 
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T. 


Pr 


n-1 


^ ( I'IP2\ n 

•^MPI ''P ■' » 


MP2 ~ -^MPl 


T 


tlPl 

n-1 

T (IWl) ^ 


^IP3 ~ MP2 


m2 


and 


T 


BIPO 


P P-.~l 

m ( IffO \ n 
■^MP3 


(3.108) 

(3.109) 

(3.110) 


Empirioal relations used to find ont the oiithalpy of 
stoan at the- exit of each M.P. turbine stage am prosonted 
below. 


%n?i 

^^EP3 


= -0.35 + 0.52 + 571.8, 

= -0.5 P^2 ^I.IE>2 5S0.4, 

= -0.636 Pj/jp^, + 0,5 + 584.5 


and 


E 


■MPO 


Efficiency correction factor foi’ each stage 
is aalculated using the folloxfing oquations ; 




\U 


'*^E1 


¥■ 


(frr~)(l - Wt.) 


>¥ 


R 


(3.111) 

( 3 . 112 ) 

(5.113) 

(3.114) 


f 


(3.115) 
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ri]yjp3 - 

and 

Actual ontlmlpios arc detominod as follows ; 
^%®1A " %!PI “ ^MPl ^^MPI “ '^%iPl^" 

^^MP2A " %P1A ” ^I-IP2^%U ^iyn?2^^ 


( 3 . 121 ) 

( 3 . 122 ) 

( 3 . 123 ) 

( 3 . 124 ) 
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^MP3A ^iyn?2A “ ^MP3 ^^M]?2A “ 

and 

^MPOA “ ^MP3A ~ ■Hf^o ^%IP3A “ 


(3.125) 

(3.126) 


Power generated "by each turbine is given by the following 
expressions 


P0¥2 = 0,0042 Wjyjpj “■ (^>.127) 

P0¥^ = 0,0042 Wj,£p2_ (3.128) 

P0¥^ = 0.0042 ''^])CP3^^]yip2A “ ^1^23^ (3.129) 

and 

PO¥j = 0,0042 "'^I>CP4^^'^MP3A “ ^%€P0^ (3.130) 


3 . 4.9 low pressure turbine 

Exit from this turbine goes to tlio condenser. Extractions 
are lumped into one extraction [ 9 ] and is assumed to take place 
at the end of this turbine. Bate of change of l.P, turbine 


exit flow is given by the expression 


d f,r \ _ (1 - ^ 3 ) - \ pp 


Actual efficiency of the l.P. turbine is 


(3.131) 
calculated as 


below 

^IPA "" ^IP %5 


(3.132) 
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WllCZ'C 


" GificiGncy correction factor is given by 


(1!^) " 


11 


R 


P5 ~ “"W "“T” 

)(1 - 


(5.135) 


ActuLil onttolpy of stean leaving the l.P. turbine and 
the power produced by the L.P. turbine ai’e given by the follow- 
ing oRuations, 


IT 

•^^•LPOA 

and 


^'^FIPOA “ ^ 


IPA “ %PO^ 


P 0 ¥g - 0.0042 ^LPA^%IP 0 A 


Hu?0^ 


(3.134) 


(3.135) 


3.4.10 Condons or [9] 

Pleat flow rate from tho notal to the cooling ii®.tor is given 
by 




(3.136) 


An on orgy balanco of tho fluid inside the condonsor is 
used -bo caloulato the enthalpy of tho condonsate exiting from 
tlie condonsor. The exprossion used is 


^^00 ^^IPO “ ^'^Sl'/^'^LPO 


(3.137) 


Energy balance on the metal wall yields the following 


Q 


SM ~ 


^'CM it ^'-^CM^ 


(3.138) 


equation 
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3,4.11 Low pressure heaters 


The nodel developed for the L.P. hoater ho. I is ropresonted 
by the following equations. 


and 


are : 


and 


"'’^Li-asM “ ^LHiMw “ ^^4 rt 

(3.139) 

^'‘lhikw ^HMva ^\ni 

(3.140) 

^%H10 ^'^CO ‘^IHim-Z'^^^LHl 

(3.141) 

The eq^iations derived to ropresont the 

« 

a 

L.P. heater ho. 

^’^LH2SM “■ ^LH2M¥ ^HM2 % d^" ^%IM2^' 

(3.142) 

‘'^LH2Pl¥ ^ ^Wl¥2 ^%H2 "■ ^HVJ2^ 

(3.143) 

^^LH20 ^ %H10 ■*■ ®“LH2M^’'%I2 

(3.144) 

L.P. heater ho. Ill io reprosontod by 

the following 


equ'xtioiis ; 
^LH3SM 


'^'LH3MW “ "IifI3 


0.6 




*^I-IM¥3 ^ ^^LH3 ^ ^ ^ 


and 


%i30 " %H20 ^“lH3ffl'/''LH3 


(3.145) 

(3.146) 

(3.147) 
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Tlic following oquations constitute the nodol for the 
L . i’ , ],.i oa t or W o , IV . 


and 


iiiid 



^^LH4MW - W 

(5.148) 

^^LH41W ^ 

^HMV/4 ''^■^1H4^ %I¥4^ 

(5.149) 

^^1H40 "" 

^^LH50 ^^lH4Mh/^'^LH4 

(5.150) 

L2 High p 

res sure lioators 


Equations 

derived to r opr o sent the ILP, 

heater No, I aro 


SlIIlMW "" SlM5 % hi 

(5.151) 

^inilMW 

^HM¥5 ^'%i5 " %Ih5^ 

(5.152) 

H - 

^'^M40 EHl 

(5.155) 


The following throo oquations ropresont tlio nodol for the 
H.P, heater II : 


and 


^“j;in2SM “ 

SfflTM'/ S'MS % df 

(5.154) 

'■^lll,[2iviw 

^HMW6 ^%M6 “ 

(5-155) 

^%H20 

%Hio %I 2 ]\w/''I-n-l 2 

(5.156) 


** 1,1 T.-^'fWR ’ 



Arc. Nfs. 
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5 • 5 Model Struc ture 

The structure of the non-linear ordinary differential 
equations is of the form 


X = f(x, y, u) 

where X is tho state vector given by 
““T 

X = ^P2M’ ^dr’ ^’do» “di’ 

^do» j'sl' ^sl' .rs2» ^'^s2» ^^o^ ^SMl' '■^'SM2; 

'"'rI' ^%101» %I2^J R2» %02’ ’%)2^ 


(3.157) 


'•^RM2’ ^'^riPl» ^MP2' ’'■'MPO' -^GM^ •^•IRIl’ 




%I5’ 


y, the auxiliary variables are given by 

Y = g(x7 U) 

U is the control vector defined by 


(3.158) 


G — (^’^0o> > ^''^AR ' load) 

Substituting for Y from equation (3»158-) in equation (a. 157) 
we can rewrite as 


X = P(X,U) 


(3.159) 
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wiioi 0 1* (101101:08 "tile now fuiicbioiici,! depondo-iiCG. 

To solve tiic above set of equation a coiaputor progranno 
was developed. Those equations were FJolvod on a digital 
Goiiputor (DEG :1090) using Runga-dCutta. intcgratio]! teclmique. 

The transiont rospoiisos of the conploto syston for stop clianges 
in control paranoters wore obtained. The results of the tests 
conducted arc discussed in Ohaptor 5 , 

3 . 6 Q onc lusipns 

In tliis chapter wo have pi’esontod the riodol developed for 
the ontiro systora which consists of forty first order difforen- 
tial oquations and one hundred and seivontoon algobro.ic equations. 
Ill tho pros Gilt modol wo have included all of the important 
conponcjiits of tho plant. But it was impossible to includo the 
dynamics of the ball-mills used for pulverising the coal duo to 
the lack of Guffioient data. FurtliGr^ tho ball-mills are adjus- 
t'od (in tho present plant sot-up) miiually from tine-to-time. 



CHAPTER 4 


results for subsystem I Am irUPORTAUT RESULTS FOR 
SUBSYSTEMS II AED III 


4 • 1 Initialisation 

A inatliGiaatical model dovoloped for the Panlci thermal 
power plant was presented in Gtopter 3. Ihe steady state 
values of the state variables are to be determined as a first 
step for any further study of the behaviour of the model. 
Obtaining the steady state operating variables for the plant 
dynamics is termed as initialisation. There are two methods 
available for initialisation. In the fir-st method the left haaid 
side of the equation (3.158) is equa.tod to aero. The new sot 
of equations are now solved using a suitable e3-imination 
technique to find out the steady state values. Second method 
is a simpler method in which opozaiting data are collected from 
the plant to use as a starting point in the integration scheme 
of the equation (3.158). Latter method is adopted in the 
present attempt. 

All the three subsystems were siniulatod to obtain steady 
state values using fourth order Runga-ICutta method for integra- 
tion. In our case the integration ^cvs terminated when the 
difference between the successive values of ■ tho state variables 
is of the order of 0.001. The steady state values obtained 
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for tho Gomplete system at a load of 100 IW is tabulated in the 
next chapter. Results of the individual tests conducted on the 
subsystem I, subsystem II [7] and subsystem III |'9] are presented 
below, 

4.2 Re su lts for Subsystems I 

Two tests have been conducted on the model. 20 percent stop 
decrease have been made on the fuel flow rate to get the tran- 
sient responses. A 20 percent step decrease in air flow from 
steo.dy state have been made to conduct the second test. Results 
obtained are discussed beloiir, 

(a) 20 percent step decrease in fuel flow rate 

The responses obtained are plotted in Rig. 3. As the fuel 
flow rate is suddenly decreased the temperatui-e of the flue gas 
also undergoes a sharp decrease initially. This is quite 
evident from equation (3.21). Aitei- the initial jump the 
temperature vory slowly decreases and assumes a new steady?' state 
value, \'Ie can note a simil.-.r trend in the case of the temperature 
(Tq_^jp) of flue gas entering the second pa,sB of the furnace. 

The eraissivity of the gas, eQ_ increases as the temperature of 
the gas decreases, in the operating region, which is the expe- 
cted trend for eniissivity. Temperature of the air leaving the 
air preheater II, T^p^^ gradually decreases as the fuel flow is 
decreased. Earlier we have noted the decrease in temperature 



of tho flue gas entering the second, pass of the furnace in 
wnich preheaters arc located. Due to tho decrease in flue gas 
temperature the heat transfer rate fro.ra the flue gas to the air 
also decreases. Hence ultimately the temperature of the a,ir in 
the air preheater II also decreases. 

(h) 20 percent step decrease in air flow rate 

Responses obtained on conducting a 20 percent step decrease 

in air flow rate from the steady state are presented in Fig. 4. 

Tornporature of tho flue gas, T suddenly increases to a high 

gw 

value and then gradually decreases to attain a steady state 
value. It is assumed that corapleto combustion takes place 
inside tho combustion cliamber. Thei’e is a momentary decrease 
in tho rate of flow of flue gas as tho air flow ro-te is decreased. 
Dug to this heat generation rate romainiiig the sane, heat 
transfer rate into the waterwall metal decreases. So that tho 
temperature of the flue gas may suddenly increase. When the 
system stabilises the high flue gas temperature gradually 
decreases to a steady state value. The flue gas entering the 
Second pass of the furnace is also having the same trend as that 
of the flUG gas in the first pass of the furnace. The sudden 
decrease and then the gradual incroase in emissivity is quite 
expected due to the sudden increase and then the gradual decrease 
in tho flue gas temperature. Air leaving the air preheaters I and 
II decreases suddenly at first and then gradually decreases to 
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attain a new steady state value. As we have noted earlier 
decrease in air flow rate decreases the flue gas flow rate. 

But flue gas flow rate is proportional to the heat transfer 
rate from the flue gas to the metal walls. This results in a 
decrease in the heat transfer rate from the metal wall to the 
air passing through the prehoaters . This ultimtely results 
in the decrease of temperatures leaving the air preheaters I 
and II. 

4 • 3 Important Results for Suhs . ystem II [ 7 ] 

Tliree tests have been conducted on the model for the 
subsystem II. The results obtained are discussed below. 

(a) 20 percent step increase in fuel flow rate 

This step change is effected by increasing the input flow 
rate,, ¥p. The responses obtained are presented in T'ig. 5. 

With the increase in fuel flow 3rato there xfill be an expected 
increase in steaming rate W^. Increase in fuel flow rate 
increases the waterwall temperature As increases the 

drum preBSu.re increases. A decreasing trend in the 

liBSS flow rates m^j_ and m^^ is noted. This decrease is due to 
the increase in friction losses in the downcomer and waterwalls. 
The decrease in m^^ is rapid as compared to m^^^. Thus bhere 
is a net increase in mass flow into the drum from dowiiconer 
wator-“Wall loop,. This causes an initial increase in the volume of 
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the druia water since tiie feed vrater flow rate and the throttle 
valve Opening are kept constant. This increase in drum water 
is called swelling effect. Due to tlio increase in fuel com- 
bustion, the temperature of the super heated steam T^^^ and 1^2 
increases. Since the outlet throttle valve opening area A is 
kept constant during this test, the additional steam produced lias 
to flow through the same area hence this causes a pressure build 
up. 

(b) 20 percent step increase in feed water flow rate 


The responses obtained in this aase are presented in Fig. 6, 
As feed water flow rate increases volume of drum w<ater incre;aj&os , 


Pressure in the drum increases gradually on account of the 
increase in drum water. The increased drum pressure propagates 
down-stream causing and Pg 2 "to increaso in primary and 
secondary superheaters respectively. On account of the influx 
of the feed water the mass flow rate of the iTater in the down- 


Goraor , 


increases initially till about 5 seconds due to the 


inorease in drum pressure, but decroasos afterwards due to tho 


increased friction losses, fess flow rato of two phs-se mixture 


m at wiaterwalls outlet decreases initially and is again due 
wo 

to the increased friction losses in waterwalls. The increased 


mass flow rate of steam through the throttle valve is due to the 
increase in pressure Pg 2 ^resulting from equation (3.62), 
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(c) 20 pcrcGnt step increase in throttle valve opening area, A. 


The responses for this tost are given in Pig. 7. With a 
step increase in throttle valve area there is an initial sudden 
increase in steam flow rate, This causes the rapid initial 

drop in tlirottle pressure P^ 2 * is explained by the fact tliat 


the friction losses increase in the superheater v^hile the drum 
pressure does not change rapidly (not shown in the Fig.). 

The rapid decrease in temperature 1^2 on account of the 
assumption made in equation (3.61) thnt is due to its dependance 
on hg 2 Oi.nd density of steam in the secondary superheater. T^^^ 


also gradually decreases as time elapses. 


4 . 4 Important Results for Subsystem III [ 9 ] 

Model responses for a load droio of 10 M/ arc discussed in 
this article. For obtaining these the total input flow to 
turbine, was decreased by a step to its new value at 80 M¥ . 

The various responses obtained are plotted and are presented in 
F ig . 8 . 

The curve for the Triflax reheater outflow rate, 

shows that goes doxm for about six seconds initially and 

iiUl 

then stabilises. This is beaausc of the large volume involved 
in the Triflex reheater it behaves just like a storage tank. 
Roheater outlet pressure, Pj^q^ also decreases as shown by the 
curve for Pj^q^ is dependant on the inlet pressure to tho 
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rcheaterp the flow and density. But the inlet pressure to the 
reheater, and the inflov/ rate, go down. So the docrease 
^ROl justifiod. As Pj^q^ goes down follows the 

same trend. The flow at the outlet of the M.P. turbine, 
deer cases initially and then stabilises to a new -value. Tbe pressure 
at the exit of the first M.P. turbine, Pi,,jp]_ P® affected by the 
inlet pressure to the turbiiio, the f3.ow passing through it and 
the density of the stean. The pressure first goes down because of 
the effect of decrease in inlet pressure to the M.P. turbine, 
but as flow also goes down and there is only a imrginal decreaso 
in the steam density, the pressure goes up as can be seen from the 
equation (3.99). The teraperature of the steam at the exit of the 
M.P. turbine follows a similar trend duo to its dopendanco on 


pressure and also due to the effect of second rcheater enthalpy 
variation. The curve for Pj,jp 3 _ shows the vciriation of the tempera- 
ture of the steam at tho exit of the M.P. turbine. The heat 
transfer rate from feed water heater metal to the water (both 


for L.P. Heater and H.P. Heater) shows a decreasing trend which 
can bo explained by the fact the.t as the load is decreased the 
quantity of steam extracted also goes do-i-m. As a result less 
steam condenses on the tubes resvilting in lesser heat transfer 
to tho feed water passing through these tubes. 



GI-LiPTER 5 


RESULTS FOR TliE COIffiniSD SYSTEM 

5 • 1 Steady State Behavi our 

All th-G three subsysteras are combined together to get 
tlio complete model of the thermal power system. The model Is 
simulated to got the steady state values of the various operating 
variables. Steady state values obtained are given in the 
Table 5.1, The units of the variables are described in the 
N omen cla tur e , 

The tests oonducted on tho model and results obtained in 
each case are discussed belovr, ■ ' 

5 » 2 Transient Behaviour for S o me To st In p uts 

5.2,1 20 percent step increase in fuel flox7 rate, ¥p 

Transient responses obtained are presented in Figs. 12 and 
13. (Capital letters are used for notations of po.ramotcrs in tho 
T'igs.), As the fuel flow rate is increased suddenly heat transfer 
rate from tho flue gas to the waterwall also increases rapidly. 
This loads to a sudden increase in the watorwall metal tempera- 
ture as shown by tho curve c of Fig. 13. But as time elapses 
more heat is being aarried by tho flow of wr^itor. Thus T^^j 
sottloa down to a slightly elevated temperature in due course of 
tine as the rate of flow of heat from the elevated furnace 



Table 5.1 


Op ora ting 

Operating 

Operating 

Operating 

Variable 

level 

Variable 

level 

~ 




dr 

335.6 

^ m 2 

31.3 

°ri 

94.8 

^ m 3 

29.1 

^dr 

139.0 . 

^MPO 

28,2 

°di 

67.4.7 

^^ROl 

97.2 


368.8 

‘Wi 

92.0 

in 

■ wo 

674.7 

^'•'gr 

97.0 

T 

■^Sl 

355 .4 

2 

90.2 

®S2 

529.8 

^'’^I.IP3 

88.3 

^sl 

140.9 

^'^I4P0 

84.2 

■^a2 

130.8 

m 

567.0 

'd 

105.5 

T 

■^ICP2 

556.0 


53.5 

1]V[P3 

542,0 

^ROl 

51.5 

^MPO 

536.0 

^RSl 

541.0 

^LPO 

45.0 

'■^’gr 

388.5 

^IIIW 

1683.0 

'■^’rS2 

647.0 

^Hir-M 

3104.0 

%02 

49.4 

'■^PIA 

155.7 

^''^R02 

97.2 

^P2A 

325.6 

^‘jyiRi 

33.2 


283.4 

%i2 

107.0 

% 

0,211£ 
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tGmpc-i’ature boooraes equal to the rats of heat carried away by 
i/lie flowing water in the watorwall. As increases more 

steam is generated in the watorwall tubes. The two phase mixture 
with higher steam content (and consequently lesser water content) 
will flow to tho drum resulting in a higher rato of circulation 
in the downcomer-waterwall loop. This is shown by the curves 
G and f (Fig. 13). As per assoimptioii, the steam outlet from the 
drum is considered to be through an imaginary throttle valve of 
constant opening area. This results in higher rate of steam f3.ow 
into the drum through the waterTCill and, consequently, higher 
drum pressure as shown in Fig. 13(a), As T^^ and are 
related by equation (3.47)j'T^^ vrill also have a rising trend-. 
Afterwards, as P^^ gradually stabilises T^^ also follows the 


trend of P 


dr* 


Trends of temperatures of stoam in primary 


superheater, T^^ and secondary superheater 1^2 Q'lKd their pressures 
Pg,j_ r and Pg 2 given by curves d,g,h and i respectively. As 
the steam generation is increased raore steam flows through the 


superheaters. But the tlirottle valve after the secondary super- 
heat or is kept at the same area of opening. Due to increased fuel 
f3.ow rate, the superheater temperatures -also increase becau.se 
the superheaters are receiving heat from the furnace. Conse- 
quently, it will add to the increase in steam pressures and 
tomperaturos (P^j^, P^g? ^sl ^s2^ of the suporhoaters . So 
there is an initial pressure building tendency and temperaturo . 
increasing tendency inside the superheaters. As 1^2 increases 
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initially the steam flow rate <^13 0 increases initially. 

All the parameters show a stabilising tr’cnd after the initial 
changG from the steady state values. The trends of the various 
parameters obtained in [7] are quite similar to the trend 
obtained in the present study, for the initial period of time, 
as shown in Fig, 5. In the pi’esent case as time elapses there 
is a stabilising trend. But in Srinivasan' s work there is no 
stabilising trend. This may be duo to the fact tlaat in the 
present model a closed loop for feed water flow path is used 
while for the one in [7] an open loop path for feed mter is 
assumed. Increase in is more rapid when compared with 


the increase in m Fig. 13 (o and f). This fact leads to a 
decrease in drum water volume initially. But after 22 seconds 
becomes greater than m^^ leading to an increase in the 
drum Wilt or volume. This increasing drum water volume phenomenon 


is called swelling in the drum. Gradually as tine elapses 
swelling dccroases. After the initial increase in m^^ and m^^^ 
both decrease due to increased friction losses. Our results 
in this case do not entirely mtch with the results in [7] and 
[1], Curve CO of Fig. 12 represents the trend of the high 
pressure heater II out-lct feed water flow rate. H.P, heater 
II out-let flow rate is assumed to be the sane as that of the 
inlet flow rate into the drum. Here we notice an initial 
decrease in we have seen earlier increases 
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initially. Both of the above facts load to a,n initial rapid 
increase in In the cases in [7] and [1] Demins 

constant during the test due to its open loop nature for feed 
water flow path. But in the present case also a 

variable duo to the closed loop nature of the feed water flow 
path. These my be the reasons for tho i-nisinatch. 


Steam leaving the H.P. turbine into the reheat er, 
is represented by the curve dd of Pig. 12. Pressure and 
t crap ere; ture of tho H.P. turbine, P^^j^ and T^p are showing an 
increasing trend initially as shown by curves k and 1, This 
is quite expected due to tho increase in VJp. Equation (3.85) 
shows tint when the density of the steam inside the Triflcx 
reheat or (which also receives heat from elevated flue gas 
temperature) increases, there is a decreasing tendency in 
This is shown by the curve aa of Pig. 12, But after five 
seconds increases and then gradually decreases and sta- 

bilises. Duo to the increased heat inflow rate into the relieator 


notal the pressure and temperaturo inside the reheator incroase 
as explained by the curves for Bp,Q-L ‘^HSl’ '^j-'ends for PpQ 2 _ 
and Tpg 2 . given by tho curves n of Pig, 15 and bb of Pig. 

12 respectively. The steam which flows out of the first a,ncl 


Second stages of the M.P. turbine, ^'^MP2 ' ^ 

similar trend as tlx'.t of and is shown in Pig. 12(ee, gg). 

Those are expected trends booause of the behaviour of 
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gradually decrease and then stabilise. Due to the behaviour 
of s beam flow out of the drun a,lso increases initially and 

then decreases. This fact leo.ds to an initial increase in 
primary superheater pressure and secondary superheater 
pressure Pq 2 * Similar trend is noted in the case of respective 
tempera tures T^^ and Tg2* both the superheaters pressures 

and tempera turos stabilise after some tine. Due to the rapid 
increase in downcomer inlot flow, m^^ increases initially. 
After the initial, increase decreases due to two reasons,* 
one to due to the gradual decrease in P^^, and the other is due 
to the increased friction losses. Due to the increase in 
watorwall outflow rate m ^ also ijicreasos. After about 20 
soGondo a little 'swelling* (^’ig. lie and f) is also observed. 

Duo to the momentary increase in 1^22 the stoo.o flow rate out of 
exit r cheat er also increases momentarily. When P^g decrease 

and stabilise around the steady state value also shows a simi- 
lar trend. Curves h and 1 (Pig. 11) depict the behaviours of 


P ,^2 respectively. 


Steam outflow rate from H.P. turblio behaves just like 


that of This is quite expected. Steam pressure P^p 

leaving the H.P. turbine also shows and expected increasing 
tendency initially. After reaching the maximum values P^p and 



doGroa.se and stabilise around 


steady 


state values. 


Curves 


GG and dd of Pig. 10 depict the behaviour of T^p and P^p 
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respectively. Due to the increased inflow rate into the 
Triflox rehoaterj density inside the Triflex roheater increases 
initially. This leads to an initial set back in the flow rate 
out of the Ti'iflex reheator, (Eqn, 3.63) > This set back is 

propagated to the f3,ow rates in all the conponents after the 
Triflex roheater. After the initial set back in it 

increases and then gradually decreases to stabilise around a 
steady state value as shown by ourve bb of Fig, 10, Curves as. 
and of Fig. 10 e.nd j of Fig, 11 show a sinilar trend as 

tluat of for %IIi2 ^''GsP'^ctively . Fuel flow 

rate is assumed to bo constant through out the test. V/e have 
noticed a net increase in all of the flow rates. Due to the 
increased flow rate more heat transfer from the metal wall to 
the working fluid takes place. This increased heat transfer 
rate loads to the docrease in the metal wall temperature. As 
the metal wall temperature decreases the heat transfer rate 
froiii the metal walls to the preheater I air and econoraiser water 
decrease. This results in a decrease of temperatures of air 
preheater I air and economiser v/ator as shown by curves 1 and k 
of Fig. 11. Following the sane argu.mont as in the previous 
article the bohiviours of Q^-, and can he explained. 

They are represented hy curves ff and gg of Fig. 11 respectively. 

5.2,3 20 percent step increase in throttle valve area, A 

The vai-ious responses obtained for the area ciiango is 
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plotted in Pig, 9. Throttle valve area A is given a 20 per- 
cont step increase to conduct the tost* This causes a sharp 
inoreaso in stoan flow rate, Due to the sudden release of 

steam out of the exit suporlioatoi' causes a drop in ozit super- 
hoator pressure 9d). The belxavioux’ of and 1^2 

are as shown by curves o and d of Fig. 9. As tine elapses both 
l/j and Pg 2 stabilise . As the pressures in the suporlioators 
docro.iBc the stoan floxir rate out of the drum increases 
initially. But as tine elapses shox^s a stabilising tendency. 
The rapid decrease in temperature of the exit reheat er steaia 
temperature 1^2 "to 'i^ho assuiaption made in the equation 

(5.61) tint is due to its dopendance on h ^2 density in the 
exit roheator. Due to the nearness of the throttle valve to 
the secondary superheaters, xdien a step change in area is given 
exit rohoater steam properties are more affected than that of 
Triflox roheator steam properties (Fig, 9 b and c). 1^2 flu- 
ctuates much more than that of Puc to the pressure 

decreaso in Triflox roheator T^^ also deer oases. Bxit 
stabilises more quickly tlmn 1 ^ 2 . Is increases heat flow 
rate into the xfaterwall does not increase because fuel floxf 
rate is kept coiistcuat during' the test. This leads to a 
docreaso in P^^, latent heat of evaporation increases with 
doGreaso in pressure P^^. This xfill again decrease P^^ duo to 
the lesser amount of steam generation. This decrease in P^^ 
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propagates down stnoam to cause a decrease in and Pg2 
(Pig, 9d), Triflox relieater stoari tenpcrature decreases 

iiiitit.illy and then increases gradually to stabilise. Since 
Wi increases the flow rate through the Triflox roheator 
and the flow rate through the exit roheator increases 

(Pig. 9g). But the heat inflow rate into thv::, notal wall does 
not increase. As increase heat inflovr rate into the steam 

from the notal wall increases. Due to this there is a decr’case 


in ].'’ohoator metal toII temperature. This fact leads to the 
deer case in initially (Pig, 9f). Because increases the 

flow rate through the H.P. heater No. II, , shows an expected 

inci'Gaso initially (Fig. 9h). stabilises as tiiee 

elapses. Fig. 9i shows the behaviour of I.P. heater No, I notal 
to water heat transfer rate, Due to tho increase in flow 

rate through the I.P, heaters also increases. After the 

initial increase like all other parameters ^Dso shows a 

stabilising trend. ¥e have seen that tho flow rate increases in 
most of the components leading to the decrease in metal wall 
temperatures of various components. The same argument holds 
good in tho case of air preheaters also. Temperature of the air 


leaving preheater I, decreases initially and then gradually 

inoreasGS to stabilise aromid a stoadjr state vo.lue. 
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5 . 3 QQoaonts 


The stahilising tendency of the paranoters under considera- 
tion can be more clearly observed and explained when the entire 
systeo in the closod loop nodel for the feed v/ater is considered 
as in tho presont case. Such trends vrere either only partially 
observable or absent from the incoupleto models studied in sub- 
systons II and III separately. 



CHAPTER 6 


OONOLUSIOHS AND SUGGESTIONS FOR FUTURE V/ORK 


^ 1 Con plus ions 

A inatlienn tica.! model for a 110 W\J thermal unit at Panki 
Thermal Power Station has been developed. For the sake of 
simplioity the entire theriml unit has been divid.ed into three 
ma 3 or subsy o t ems . 

i) Subsystem I consists of eoonomiBor, air preheaters, 

fuel feeding system including pulverisers cmd inrued draft 
fan, and the furnace including the burners 

ii) Subsystem 11 consists of boilei" drum^ downooxaer* water- 
walls, and primary and seconda.ry superheaters. 

iii) Subsystem III consists of high pressure, medium presmiro 
and low pressure turbinos, steam regulating system, 
throttle valve, relieators condenser and feed water heaters 

Subsystem I has been modelled in this thesis whereas the 
models for subsystem II and III have been referred from [7] and 
[9], The thx-ee systems have then been combined including 
a closed loop feed w'ator path. Thus a complete model for the 
entire thorml side of the power plant has been obtained. 
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llie entire mathomatical model Ins forty ordinary non— 
linear differential equations and' one hundred and seventeen ali-ge— 
braic equations involving non-linoar relationships between the 
parameters of the system. Some of the important input parameters 
whoso variations studied are fuel flow rate, feed water flow rate, 
and throttle valve opening area. The effect of these variations 
on the transient belxiviour of some important output parameters 
such as drum pressure, superheater pressures, steam flow rate 
through throttle valve, steam flow rate through the turbine stages 
etc. have been obtained. 

In this study a stabilising tendency of most of the sye-tam 
paramotors is observed. This can bo explained when the entire 
system in the closed loop model for the food water is oonsidered. 
Such trends wore either only partia,lly observable or absent from 
the incomplete models studied previously in the subsystems II 
and III separately. I'rom the plots, the system parametoar varia- 
tions on the boiler and turbine side can be simultaneously 
studied , 

S • 2 Sug ge stion for lut-iir o Wo rk 

The above study was carried out under open loop control 
mode and t'he automatic control schemes for variations in power 
level and. turbine speed have not been taken into account. 

The control sohomes and hoirarcl];^'- should also be considered as a 
next stop towards completing the model. The model on the 
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olGcbricsal genera 1:1011 side iirhicli is c^Hipled. Iio tlae turbine, 
should also be incorporated to study the eiitiro system. 

Obviously this is going to be more exhaustive and difficult. 

An attempt was made in riirchp 1978 to procure transient 
data by conducting 10 MV/ and 20 111 step power drops at Panki 
Thorriial Power Stationi . In these tests open loop- -c-oixtrcil mode 
was used to simulate our model. Unfortunately due to in- 
adequate and sonewliat crude instrumentation (which may be 
suffioiont for the effective operation of the plant) only 
qualitative trends for some of the parameters could ba obtained. 
Some of these trends match qualitatively with our theoretical 
results. It is hoped that tfes mox'e aocuraie study oouid be done 
by conducting actual tests on the plant with someiihat more 
sophisticated instrumentation. 

The integrated model for the entire unit Including theriml^ 

systems 

generation, distribution and control/can be simulated on a hybrid 
facility. This can bo used for generating theoretical data 
for the parameters of ohoice which can be used for future desigr 
The facility can also bo used as a trainer for the technical 
staff dirootly rosponsible for running the power station. 
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APPENDIX I 


The important specifications of the 110 IW Panki Thermal 
Power Plant 


FJanuf a c tur er s BIIEL , BID lA 

Generating Capacity 110 IW 


Steam Flow rate 

Live Steam temperature 

Pressure at the final 
superheater ou.tlet 

Fuel 

Calorific Value of the Fuel 
Ash content of the fuel 


375 Tons/hr 
540°0 

139 Ifg/cm^ 

Pulverised Coal 
4727 X Cals /Kg 
32 percent hy vr eight 


PLANT DATA 


Heating surface area 
of economiser 

2880 rs? 

Heating surface area of 
air preheater one and two 
together 

18720 

Metal weight of economiser 

62570.15 ICgs 

Metal weight of air 
preheater one and two 
together 

78620.36 ligs 

Water content in the 
economiser 

22 nP 
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MODEL CONSTANTS 


^PIOM 

^PIMA 

^PIA 

^PIM 

^P2GM 

‘^PPiyiA 

°P2A 

Cp2iv[ 

GOBI'' 

^EM¥ 

^ims 

'■^002 

^N 

^HgO 

^throttlG 

CONST 

0 


0.8099 

0.8022 

0.25 

0.164 

1.055 

1.1627 

0.25 

0.164 

1.063x10 

2361,36 

0,6125 

1.01 

0.16608 

1.0 

0.878 

0.5899 

0.04919 
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NOMENCLATURE 



a 


K 

AUS 

Ai jAg 


GOA 

Cp 

G 

im 

^S2 

CONST ; 

0 ^ • a e ft 0 ^ o 


Cross sectional area of down-ooiaer , (nr^ ) 

Absorptivity of fieJB gas 

Plow area available, for a partially 
opened valve (m^) 

Cross sectional area of waterwall tube&(na^ ) 

Watorwall effective surface area,{in^) 

Equivalent areas for I and II reheater 
tubes respectively (m^ ) 

Constants used in flow» pres-suro relations 
for governing sttige 

Constants used in the equations for flow 
Area of the throttle valvoj (n'^) 

,) 

Specific heat of fi.a© gasas^ -(kco-ls/kg) 

Constant for pipe flow from drum to 
superheater 

Calorific value of pulverised coal#-(koals/lcg) 

Coefficient of heat transfer from super 
heater metal to steam in primary and 
secondary superheaters. 

Coefficient of heat transfer between 
primary superheater metal and furnace 

Coefficiontefheat transfer between socon- 
dary superheater metal and furnace. 

Specific heat of waterwall^^and secon- 
dary suporhoater metals respectively 
(kcals/lcg °0) 

Cooff iciontsf’heat transfer between 
waterwall metal and two phase mixture 

Model constants 



Tfiii 


°EGM 

^P1M'^P2M- 

^EM 

^P1A»^P2A ° 
C 

EW 

^PIMA’ 

^P2AMA 

^TMl’ 

^TM2 

^RMSl’ 

^RMS2 

C ; 

cw 

‘^HMW1'^HMW2 
^HMW5 '^my[w6 

.6E ; 

'^G'®C02’ “ 
%0 

f 


CoGfficients of heat transfer h.e-trvreen 
flue gas to metal of air preheater I, 
air preheater II and economiser 

Specific heats- of airpreheater I and II 
and economiser (Kcals/K^ °C) 


Specific heats of Air inside airpreheater 
I and II ( Reals /IQg OQ) 

Specific heats of water inside the 
economiser (Iicals/li;g ^0) 


Coeff icientc^heat transfer hetwoen air- 
pr'sheater I and II and air 

Coef f icionte{-hoat transfer between econo- 
miser metal and water 

Specific heat of the metal in reheater I 
and II (Reals /Kg ^0) 

Constants used in the model 


Constant 


Constants used in the heat flow equations 
in the four L.P. heaters rojspectively 


Constants used in the H.P heater model 

Piamete-r of Downcomer and waterwall tubes 
respectiTelyo (n) 

Correction factor used in emissivity 
equation 

Emissivity of f-i-^ii gas, carbon dioxide 
and water vapour In the combustion 
chamber respectively 

Vector function 



ix 


f f 
1? TT 

^Rl’-^R2 




%S’^'^HP 


^'’i'^’^di’ 

^s2’^^sTl, 
h, 


'ST2 


%01'%02^ 

%IR?/‘^MP1' 

^■'^MPO 

%D?1A,%"2A * 
^%[P3A’%^OA 
%PO»%POA ’ 


%jH10'%i20 • 

% H 30’%140 

%nil0'%IH20 

J 


j Ic 2 


Pr let ion factors in tli-,: downcomer, water- 
\'ra.ll, Triflex and exit rolieator respectively 

2 

Acceleration due to gravity (n/sec ) 

Enthalpy of '^eedwater, att c ipex’ation spray, water 
in the drum, steam in the drum, steam leaving the 
drum, water at downcomer inlet, steam at the 
outlet of water walls, steam at the exit of thrO“ 
ttle valves and steam at the exit of governing 
stage. (Kcals./Eg) 

Enthalpies of steam leaving drum, water at 
downcomer inlet, steam in primary super- 
heater, steam in secondary super-heater, 
steam at steady state operating temperature in 
super-hea t er s ( Zeals . /Zg ) 


Enthalpies of rehea,ter steam (ideal and 
actual), steam at the exit of triflex and exit 
reheater, att empeivition water, ideal enthalpy 
of steam at the exit of each M.P. turbine, 
(Kcals./Eg) 


Actual enthalpy of steam at the exit 

of each M.P. turbine (Zeals . /Ilg) 

Ideal and actual enthalpy at the exit of 
L.P. turbine (Zeals. /Eg) 

Enthalpy of water at the exit of four 
l.P. heaters and at the exit of two 
H.P. heaters. (Zeals. /JIg) 


Conversion factor for converting work to 
heat energy 

Callender's empirical constants. 



X 


K jKo 

1^2 

K3.K4.K5 




Constants used in the nod el 


Vortical length of downconer, flow length 
of downconer j, flow length of water wall (m) 


^1 '^2 
Mf j M2 > 

^P2M’'^EM 


Equivalent length of tubes in the triflex 
and exit rcdaeater respectively (m) 

Maos of primary sup.hoator notal, 
secondary sup. heater netal reheaters one 
and two. Air preheaters onca and two and 
ocononiser netal. (Kg.) 


n „ , n . , ; 

00 ’ ri ’ 

n 1 . s n , , n , 
di’ do’ wo’ 


M 


Rl 


M. 


R2 


^^P1A’’^P2A’ 


M 


E¥ 


Hass flow rate of feed water, stean going 
out of drun, water at dox'jn coner inlet, 
water at downconer outlet, two phase 
mixture at waterwall outlet, effective 
stooii nnss at reheaters one and two 
respectively. (Kg/sec.) 

Rcsidont mass cf air in airprohoaters 

one and two and water in the e-concmisor (Kg.) 


^SS ” 

Stoaan pressure beforo throttle valves (Kg/m ) 

^dr’^do’ “ 

P , P 
■^al'-^a2 

Pressure at drum, downconer outlet, 
pressure in the prina.ry and secondary 
sup„ heaters respectively (Kg'/m2) 

Pg.Hl(I=1.2,3) 

s Nozzle ring pressure in the three streams (Kg/: 

%P ° 

Steam prcs,s^vre at the exit of curtis 
stage (Kg/n^) 

^CR 

2 

Gold reheat steam, pressure (Kg/n ) 

^"r 01’^R12’ ” 
\pi’%lPl, 

%P2 ’ %P3 ^ 

%P0’^LP0 

Stoaxi pressoire at the exit of triflex 
and exit rehoater. Inlet steam pressure to 
M.P. turbine, pressure at the exit of 

4 M.P. turbines and condenser pressure 
respectively. (Kg/m^) 



P0¥2jP 0¥^, ; Power prorluced by each of the M.P, turbines 
POW^ ,P0¥-^ respectively (M¥) 


POV/ 


6 


Power produced by L.P turbine (IW) 


%\>1 


Heat flow rate to stea,n in the priaary 
and secondary supheater, two-“pho,se nisture 
in waterwalls, primary .sup. heater netal 
from fila.^ gases, secondary sup. heater 
netal and waterwall cetal from gases 
respectively (Kcals/sec.J 


^Bl’^S2 


Net heat content of stean in the prinary 
and secondary suporheaternrespectively .CKoalS.) 


^RS1’\S2’ 

'^M¥ 


Heat flow to triflex and exit reheater 
netal walls, heat flow to trifle*? and exit 
reheater stean and heat flow to the 
condenser cooling water, (Kcals/seo. ) 


^LHIMW ’ 
^LH2MW’ 

\h3M¥’ 

^LH4MW ’ 

^HH1M¥" 

^HH2M¥ 


Heat added to feed water in the four I.P. 
heators and two H. P. heaters (Ecals'/seo. ) 


QpiQM»<^P2GM ° 
^EGM 


Heat flow to netal of Airpreheater I and II 
and econoniser netal fron gas. (Kcals/sec.) 


^PIMA’^^PPMA 

'^BMW 


Heat flow to the air in Air preheater I 
and II fron the netal and heat flow to 
the econoniser water fron the econoniser 
netal (Kcals/sec.) 


Stefan Boltznan constant u,sed in radiative 
heat transfer calculations. 



zii 




Tine- constT,nts of H.P. turbine, four 
M.P. turbines, L.Po tu.rbine unci “triflex 
reheuter und Exit rohouter 



Tine constant in the sup. heater 


'^gw 5312 


^R1’^R2' 


Saturation tonp.in the drun, Terpcrature 
of stean in the prinary a,nd secon''''ary 
sup. heater, notal of prinary and 
secondary sup. heater, f^oB gas in the 
furnace, steady state operating tenp. 
of stc-an in the primry and secon^^ary 
sup. heaters (^C!) 


Constants used in tho nod el 


^RMl’'^EM2’\si^ 

^RS2’^]y[Pl’^MPl' 

^MP2’^MP5’^MP0 


Tenp. of stean before throttle valve, 
stean at the exit of curtis stage and 
cold reheat stean. (<^G) 

Tenperataire of tube netal of trifler 
and exit reheater. Stean in trifled and 
exit reheater, stean at inlet to 
M.P. turbine, stean at the exit of 
four M.P. turbines respectively. (°G) 


^iiro'^CM''^cw’ 

m rn m 

■^HMl ' HM2 ' •'HM3 ’ 

^H\a’^HW2’^HW3’ 

^HW4’^FIW5*^m-/6 


Tenporature of condenser tube netal 
and cooling water, tube netal of 
four l.P. heaters and two H.P. heaters. 



Tonp. feed water in four L.P. heaters 
and two H.P. heaters (°G) 


^PIM'^PIG'^PIA’ 

^P2M''^P2G'^P2A 

rn m rn 


Tenperature of Air prohe-ater I netal, 
fiift^ gas and air, Air iDreheater II 
netal, fitaife gas and air, Bconoriser 
netal, f-^afe gas and wa.ter. (°G) 



ziii 


V 

V 

V 


s2»\M’'^drw» 

drs 


^Rl’\2 


1 ' o 




'G, 


^^ASH 

^DS1’^DS2 


(¥n,) ,(¥„) , 

^ V ^ N 
^GR’ ^R/\ 


¥. 


Volune downcoLier tubes, water wall 
tubes, drum, prinary and secondary 
superheater, water wall tube netal, 
drur water, ste.a^ in the drun 
res]pectivGly, 

Effective volume of steam in Trifler 
and exit reheat ers (ci^) 

Mass flow rate of steam at secondary 
and primary sup. heater outlet. (Rg/sec) 

Eiiaii flow rate aud Air flow rate and 
fila^ gas flow rate (Kg/sec.) 

Ash content of pulverised coal.C/) 

Mass flow rate of attemporation spray 
I an^ II (Kg/SGC.) 

Flow rate, through the throttle valve, 
through the nozzle ring, flow rate 
of cold reheo.t steam and ratrjd flow 
of the turbine. (Kg/seo.) 

Ratio 'f ¥jj n.nd 


^RI2’^R01’ 

^^R02’\tMR' 

¥ 

^MFI 


Plow rate of steam air exit reheater 
inlet, at the exit of trifleir and exit 
reheat or, attemperation added to 
steam before M.P. turbine inlet, 
steam, entering M.P. turbine. (Eg/sec. } 


^M.P1’^*^MP2’ 
^Mr5 ^ ■'^MPO ’ 

^''^C¥ 


Steam flow at the exit of the four 
M.P. turbines respectively,, Plow to 
condenser, quantity of cooling 
water in the condenser. (Eg/seo.) 


l’\ 2’ 
^%H1’%I2 


Quantity of feed water passing through 
the four L.P. heaters a/nd two H. P. 
heaters. (Eg/soc.) 



xi,r 




WO 


drs^ /*wos^ 

di ’ i^do ' f\ 

? si»-Ps2»fwM» 

STi'J’sa^ i^air 


Density of water in the rlrua, stean in 
the clrun, wo,ter at downconer inlet 
and outlet two phase nixture at water- 
wall outlet, stean in the primary and 
secon'-hary sup. heaters, waterwall 
netal, stea-n at stoa^dy state operating 
temperature in sup. heater I and II and 
air from the atmosphere respectively. 


Jss* ?THI’ Jhp' 

J®R1’ pR2'i^MP2, 
f MP2, j* MP3 


Density of stean before throttle- valve 
stean. in the nozzle ring, stean a,t 
the exit of governing stage, stean in 
the trifler and exit reheater, stean j 

in the three M.P. turbines respectively{Kg/ni ) 


D. » DhP’ %P’ 
T)p2' ^I'3' 


Burning efficiency for combustion, 
efficiencies of H.P., M.P. and L P. 
turbines, efficiency correction 
factor, for H.P. turbine, for each of 
the four M.P. turbines and correction 
factor for l.P. turbine 


^PA’ \vV 
^MP2’ ^Mr3’ 
'’^.MP4’ '^LPA 


Actual efficiency of H.P. turbine, 
Actual efficiencies of each of the 
four M.P. turbines an'"'' actual 
efficiency of L.P. turbine. 
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